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ABSTRACT 

Context. This paper is part of a series involving the AMIGA project (Analysis of the Interstellar Medium of Isolated GAlaxies), which 
identifies and studies a statistically significant sample of the most isolated galaxies in the northern sky. 

Aims. We present a catalogue of nuclear activity, traced by optical emission lines, in a well-defined sample of the most isolated 
galaxies in the local Universe, which will be used as a basis for studying the effect of the environment on nuclear activity. 
Methods. We obtained spectral data from the 6th Data Release of the Sloan Digital Sky Survey, which were inspected in a semi- 
automatic way. We subtracted the underlying stellar populations from the spectra (using the software Starlight) and modelled the 
nuclear emission features. Standard emission-line diagnostics diagrams were applied, using a new classification scheme that takes 
into account censored data, to classify the type of nuclear emission. 

Results. We provide a final catalogue of spectroscopic data, stellar populations, emission lines and classification of optical nuclear 
activity for AMIGA galaxies. The prevalence of optical active galactic nuclei (AGN) in AMIGA galaxies is 20.4%, or 36.7% including 
transition objects. The fraction of AGN increases steeply towards earlier morphological types and higher luminosities. We compare 
these results with a matched analysis of galaxies in isolated denser environments (Hickson Compact Groups). After correcting for the 
effects of the morphology and luminosity, we find that there is no evidence for a difference in the prevalence of AGN between isolated 
and compact group galaxies, and we discuss the implications of this result. 

Conclusions. We find that a major interaction is not a necessary condition for the triggering of optical AGN. 
Key words, galaxies: evolution - galaxies: interaction - galaxies: active - surveys 



1. Introduction 

It is widely accepted that galaxy evolution is strongly influenced, 



or even driven, by environment (nurture; see e.g. Park & Choi 
Liu et aL]|2012 and references therein). Galaxy interac 



2009 



tions are thought to induce nuclear activity by removing angu 
lar momentum from the gas and, in this way, feeding the cen 
tral black hole (|Shlosman et al.|1990}|Barnes & Hernquist|1991 



|Haan et al. 2009| |Liu et al.|201 1) >. Hence, a higher rate of active 
galactic nuclei (AGN) is expected among interacting galaxies. 

However, different studies have yielded contradictory re- 
sults: some studies find a higher local density of companions 
near galaxies hosting an AGN or a higher prevalence of AGN 
in interacting galaxies (|Petrosian|l 982; Dahari 1985; MacKenty 



1990; Rafanell !~eTal.||i995| |19971 |Alonso et al.||20071 |Ellison 



et al.|20lTj|Liu et al.|2012|l, w hile others find no excess, or only 



1995 



a marginal excess (Bushouse 1986; Laurikainen & Salo 
Schm itt||200T] |Miller et al.||2003t |Ellison et al.||20081 |Li et al'.| 
2008). Schawinski et al . I d 2009 1120 1 01) found a clear relation be 



tween some AGN and the signatures of a previous merger. In re- 
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* Full tables 1 7, A.l A.3 are only available in electronic form 
at http://amiga.iaa.es/ and at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
|bin/qcat?J/A+A/545/A15| 



cent studies the relation between X-Ray AGN and environment 
was also explored (Silverma n et al.|2009 Haggard et al.| |2010 



Tass e et al.|201 1) , finding that this relation may depend on the 
mass of the host galaxy ([Silverman et al. 2009 ). Another possible 
source of discrepancy involves the methodology used to estab- 
lish the presence of an AGN. Subtraction or non-subtraction of 
the stellar component from the nuclear light may lead to differ- 
ences in the nuclear emission classification. This is particularly 
true for weak emission lines, which might be affected not just 
by a low signal-to-noise ratio, but also by details of the stellar 
spectrum (Ho et al. 1997| l. The use of different classification cri- 
teria to evaluate whether the nuclear emission is powered by star 
formation (SF) or by an AGN also complicates the direct com- 
parison of results found in different studies. 



The AMIGA project ( |Verdes-M onteneg ro~et al.| |2005| 
Analysis of the interstellar Medium in Isolated GAlaxies; 
http://amiga.iaa.es/) was initiated to separate the effect of galaxy 
interactions from the intrinsic evolution in a galaxy. AMIGA 
provides a panchromatic characterisation of a well-defined and 
statistically significant sample of isolated galaxies. AMIGA can 
be described as a vetted or value-added catalogue, based upon 
the Catalogue of Isolated Galaxies (CIG; Karachents eva| 1 973) , 
which is composed of 1050 galaxies compiled using an isolation 
criterion that implies that the galaxies have probably been unper- 
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turbed for » 3 Ga ( Verdes-Montenegro et al. 2005 £ Dwarf com- 
panions (B>-18 to -15 depending on distance) are not (cannot 
be) excluded, but the degree of isolation from major companions 
(greater than 10% of the mass of the primary) was re-evaluated 
and quantified for each galaxy in terms of both the local number 
density of neighbours and tidal strength ( Verley et al.|2007b|a) >. 

Multiwavelength studies of the AMIGA sample (data re- 
leased via a Virtual Observatory interface) have placed special 
emphasis on structural parameters as well as different compo- 
nents/phases of the interstellar medium. Our studies suggest that 
the most isolated galaxies show different properties than less 
isolated field samples. AMIGA early-type galaxies are usually 
fainter than late-types in the B-band and most spirals in our 
sample appear to host pseudo-bulges and not classical bulges 
([Verdes-Montenegro et aLl|200"5| |Sulentic et aTJ|20"06l [Durbalal 



et al. 2008b i. AMIGA spiral galaxies are redder than similar type 



galaxies in close pairs, showing a Gaussian distribution of the 
(g-r) colours with a smaller median absolute deviation (almost 
half) compared to galaxies in wide and close pairs. Properties 
such as far-infrared (FIR) and radio continuum emission, which 
usually show enhancement in interacting galaxy samples ( |Darg| 
et al.||2010| ), show levels in AMIGA that are at or below those 
measured in any other galaxy samples: This has been shown for, 
Lfir ( |Lisenfeld et al.||2007) , L lt4GHz ( |Leon et al.||2008), radio- 
excess above the radio-FIR correlation (0%; Sabater et al. 2008 ; 
Sabater 2009), the Hi asymmetry (|Espada et al.|201 1 1, and the 



molecular gas content ( |Lisenfeld et al.|201 AMIGA galaxies 



are also located in the lowest environmental densities (Verley 



et al. 2007a i and do not present any morphological signatures of 
interactions ( |Sulentic et al.|2006| >. 

Here, we present a catalogue of nuclear properties for 
AMIGA galaxies obtained from optical spectra provided by the 
Sloan Digital Sky Sur vey (SDSS; |York et al-pOOO} |Adelman- 



|McCarthy et al.||2008| l in its 6fh Data Release. In Sect. |2j we 
present the sample of galaxies used in this study and describe 
the compilation and reduction of the data, including the subtrac- 
tion of the stellar contribution from the nuclear spectra and the 
measurement of emission lines. The method used for the classi- 
fication of the nuclei, which takes into account the presence of 
censored data, is presented in Sect. [3] together with the final clas- 
sification of the nuclei and its relation with the properties of the 
galaxies. A comparison with galaxies in a denser environment 
(compact groups) is analysed in Sect. |4] Finally, we present a 
summary of the work, followed by a discussion about the effect 
of major- mergers in the triggering of nuclear activity and our 
conclusions, in Sect. [5] 



2. Data 

2.1. The sample 

We adopted the AMIGA sample of isolated galaxies as the 
nurture-free base sample for our study. The isolation criteria for 
AMIGA galaxies imply nearest-neighbour crossing times of » 3 
Ga or more ( Verdes-Montenegro et al.|2005| |Verley et al.|20 07a; 
Argudo-Fernandez et al. 2012J in prep.). This, and the proper- 
ties of the sample presented in the introduction, lead us to hy- 
pothesize that AMIGA galaxies were not affected by major tidal 
interactions during the last part of their lifetime. 



1 We will follow the recommendations for units of the IAU Style 
Manual (Wilkins 1995). Hence, we use the term annus, abbreviated 
as "a", for year (a - annus - year; Ma - Megaannus - Megayear; Ga - 
Gigaannus - Gigayear). 
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Fig. 1. Result of a (V/V m ) completeness test for the CIG SDSS 
sample (n = 353, red triangles) and the CIG full sample (n = 
1050, blue circles). 



The SDSS provides a large homogeneous database of optical 
spectra from which line ratios can be extracted and quantified 
to investigate nuclear activity. We used Data Release 6 (DR6), 
which was the one available when we started the study. We 
found a total of 549 AMIGA gala xies in the SDSS-DR6 database 
( Adelman-McCarfhy et al. 2008 1 with nuclear spectra for 362 of 



the galaxies. We found that spectra for 353 galaxies could be 
used for the present study as explained in Sect. |2.3| Henceforth, 
we refer to this subset of n — 353 isolated galaxies as the SDSS 
sample. 

We estimated the completeness of the SDSS sample using 
a (V/V m ) test (as in |Verdes-M onteneg ro et al.|2005| > and found 
it to be even more complete than the full AMIGA sample up 
to a limiting magnitude of mg = 15.0 (see Fig. [TJ. This re- 
flects the completeness and uniformity of data within the region 
of sky covered by SDSS. For studies that require statistical sig- 
nificance and completeness we selected those galaxies with nig 
in the range 11.0 to 15.0 ( |Lisenfeld et al.||2007 1. We excluded 
galaxies flagged as interacting in Sulentic et al. ( |2006 marked 
with the parameter I/A = yes) as well as local dwarf galaxies 
CIG 663 (UMi dwarf) and CIG 802 (Draco dwarf). We refer to 
this sample, composed of n = 226 galaxies, as the SDSS com- 
plete sample. 



2.2. Morphology and luminosity 

To study the relation of the nuclear activity type with the proper- 
ties of the host galaxies we also used morphological data, optical 
luminosity and near-infrared (NIR) luminosity. The morphologi- 
cal data are based on work presented in |Sulentic et al.| ( [2006[ ) and 
Fernandez Lorenzo et al. (20 L2]l. Optica l luminosities (Lb) were 



presented in Verdes-Montenegro et al. ( 2005| l and were refined 
based on new distances for AMIGA galaxies and a careful revi- 
sion of the velocities and corrections (Fernandez Lorenzo et al.l 



|2012 i. The infrared luminosity was derived from the K s -band 
(2.159 //m) magnitude obtained from the Two Micron All Sky 
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Survey (2MASS) extended source catalogue ( |Jarrettetal.|2 000). 
The Lk was obtained from the total (extrapolated) K s flux con- 
sidering that Lk = 47rd 2 v/k,(v), where d the distance to the 
galaxy, v the central frequency of the K s band and fy^iy) the flux 
density derived from the magnitude in the catalogue. Finally, the 
luminosity was normalised by the solar luminosity in K s band 
(L K 0). After a visual inspection we found some galaxies as- 
sociated with two 2MASS sources within our search radius. In 
each case, the flux of the brighter one is the one correspond- 
ing to the whole galaxy, hence it was selected for these galax- 
ies ( |Lisenfeld et al.|201 l| l. Although a mass-to-luminosity ratio 
could be applied to derive a total stellar mass ( |Cole et al.|2001| i, 
it depends on the star formation history of the galaxy and on the 
initial mass function considered ( |Bell & de Jong 2001 1 Bell et al. 
2003) >. Hence, we directly used Lk as an estimator of the stellar 
mass content of the galaxy. The optical and infrared luminosi- 
ties, as well as the morphological classification, are presented in 
Table Q] 

Additional data on the nuclear activity of AMIGA galaxies 
were presented in Sabater et al. (2008 ) and Sabater (2009). The 
classifications in those papers were taken from the literature, and 
obtained from the study of the radio-FIR correlation and the FIR 
colour. 

We tested whether the SDSS complete sample represents 
well some fundamental properties of its parent sample (the 
AMIGA complete sample, |Lisenfeld et al.|2007[ ). In Fig. [2] the 
cumulative distributions of distances and luminosities are shown 
(left and right panels, respectively) for the AMIGA complete 
sample and the SDSS complete sample. To quantify the com- 
parisons we used a Kolmogorov-Smirnov two-sample test (K- 
S ; |Kolmogorov| 1 933 Smirnov|1936[ ) to estimate the likelihood 
that the SDSS and AMIGA samples are derived from the same 
(undefined) parent population. The K-S test significance levels 
are 0.15 and 0.69 for comparisons of distance and source lumi- 
nosity, respectively. These values are sufficiently high to support 
the hypothesis that the two samples come from the same parent 
population at a 10% significance level (a = 0.1, which is widely 
used in the literature). We also applied a^f 2 test to compare the 
morphology distribution of the two samples. This test is espe- 
cially suitable in this case because of the discrete nature of the 
morphological classification. Fig. [3] presents the morphological 
distribution of classifications. As we can see in the figure, and 
as the x 1 test confirms, the distributions are very similar in all 
cases. The significance level is higher than a = 0.10 (p-value 
= 0.484). Therefore, we have no reason to reject the hypothesis 
that the SDSS sample follows the same distribution as the parent 
sample at a significance level of 10%. We therefore conclude that 
no clear bias exists in the SDSS complete sample relative to the 
AMIGA complete sample and we consider the SDSS complete 
sample used in this study to be representative of the AMIGA 
complete sample. 

2.3. SDSS spectroscopic data 

Spectra were found for n - 398 AMIGA galaxies in the SDSS 
DR6 catalogue. Nuclear spectra were available for 362 of these 
galaxies and non-nuclear spectra were discarded for our study. 
We downloaded the corresponding spectra and compiled the flux 
density, its error and the error mask for each wavelength. 

We extracted measures using an alternative reduc- 
tion pipeline of the SDSS data (called specBS). This 
Princeton University product is publicly available at 
http://spectro.princeton.edu/ where the spectral parameters 
are determined via a comparison with carefully selected 
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Fig. 2. Cumulative distributions of distance (left) and optical lu- 
minosity (right) for the AMIGA (solid line) and SDSS complete 
samples (dashed line) used in the Kolmogorov-Smirnov two- 
sample test. 
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Fig. 3. Binned distributions of morphological types for the ob- 
served (open bars) frequencies for the SDSS complete sample 
and scaled (grey-filled bars) frequencies derived from the CIG 
complete sample. These frequencies were used to perform a^ 2 
test. 



templates of stars, galaxies and quasars from the SDSS and 
ELODIE surveys (Moultaka et al. 2004 1. This resulted in im- 
proved measures of the velocity dispersion and the redshift. The 
velocity dispersion could be reliably determined for 42% of the 
galaxies in our sample using the standard pipeline and 98% of 
the galaxies using the Princeton pipeline. Redshift discrepancies 
between the two procedures were negligible (there are only 
three galaxies with a difference of more than 100 km/s) although 
the uncertainty in the redshift measurements is lower using 
the Princeton pipeline. As a final test, we visually checked 
all spectra and found that some of them had regions in which 
spectral data was missing. In some of the cases the spectral gaps 
did not affect the regions of interest for this study. In only one 
case (CIG 347) the entire spectral data were corrupted. 

After flagging galaxies with bad spectral regions, we have 
n = 353 galaxies for our study of nuclear spectra in the AMIGA 
sample. The spectroscopic data are presented in Table [2] In the 
column "Comments" we marked galaxies with bad spectral data. 
These nine flagged galaxies were not included in our study to 
avoid inhomogeneities in the data sample. K-S two sample tests 
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Table 1. Near-infrared and optical luminosities and morphologies of CIG galaxies. 1 



fir 1 


f 2 


l°g(Mc) 


code"* 


lOg(-LB) 


morpho. 4 




Jy 


[Lk.q] 




[Lq] 


RC3 


l 


0.0548 ± 0.0020 


11.227 ±0.016 


i 


10.57 


5.0 ± 1.5 


2 


0.0087 ± 0.0010 


10.407 ± 0.050 


i 


9.84 


6.0 ± 1.5 


3 


0.00539 ± 0.00053 




-i 




4.0 ± 1.5 


4 


0.2855 + 0.0021 


10.9777 ± 0.0032 


i 


10.28 


3.0 ± 1.5 


5 


0.01320 ±0.00043 


10.686 ±0.014 


i 


9.95 


0.0 ± 1.5 


6 


0.0129 ±0.0013 


10.205 ± 0.042 


i 


9.80 


7.0 ± 1.5 


7 


0.0197 ± 0.0012 


11.270 ±0.027 


i 


10.35 


4.0 ± 1.5 



1 AMIGA sample (n = 1050). 

2 Units: log(L K ) given in units of K s -band solar luminosity (L k q = 5.0735 x 10 25 W); log(L B ) given in units of solar bolometric luminosity 
(L Q = 3.842 x 10 26 W); Uy = lO^Wm^Hz-'. 

3 Code for the infrared luminosity: -1 if the galaxy does not have distance data; if the galaxy does not have NIR data; 1 if the galaxy has 
NIR data and distance; 2 if there are 2 NIR sources for the galaxy, the brightest one is selected (see text). 

4 log(LB) and morphology data from Fernandez Lorenzo et al. 1 2012 1. 



were performed among the raw sample, the clean sample and 
the flagged galaxies for the distance, optical luminosity and B 
magnitude. Throughout, the p-values were higher than 0.8. A 
X 2 was applied for the morphology and a p-value of 0.99 was 
found. The properties (luminosity, morphology and distance) of 
the flagged galaxies are comparable to those of the remaining 
galaxies; consequently, no bias is expected to be introduced by 
removing them. 

2.4. Stellar populations 

The spectrum of a galaxy includes contributions from both stel- 
lar (absorption) and gaseous (emission) line components. Both 
are of potential interest but they must be separated to draw any 
quantitative conclusions. The stellar absorption spectrum must 
be extracted to perform a reliable fitting of the emission lines 
(Ho et al. 1997 ). Most methods used to determine the stellar pop- 



the average of Z and log(age) at the level of up to ~ 0. 1 - 0.2 dex 
Cid Fernandes et al.|2005||Cid Fernandes & Gonzalez Del gado 



ulations are based on templates of stellar absorption-line spectra 
that are used as the basis to fit the stellar component (e.g. |Ho 



etal. 1997 Engelbracht et al. 1998 ; Kauffmann et al. 2003b Cid 



Fernandes et al.|2005||Hao et al.|20b5) >. The library of template s 
should be large to allow the fitting of stellar templates that re- 
flect a wide dispersion in galaxy properties such as metallicities, 
ages and velocity dispersions. From among the available codes 
we selected Starlight ( |Cid Fernandes et al.||2005| l, not only be- 
cause of its ability to accurately determine the properties of the 
underlying populations, but also because it is publicly available, 
easy to use, and well documented. 

Starlight is designed to fit an observed spectrum with a 
mode l ( |Cid" Fernandes et al.|2005HMateus et al.|2006HAsari et al.| 
|2007| i composed of spectral components from a pre-defined set 
of input spectra. The spectral base can be made up of observed 



template spectra, evolutionary synthesis models ( Cid Fernandes 
|et al.||2004[ |2005|>, individual stars, e tc. We used 45 synthetic 
spectra from Bruzual & Chariot (2003) as the spectral base, with 
three metallicities (Z) and 15 different ages. The 3 metallicities 
are: 0.004 (0.2 Z ), 0.02 (Z ) and 0.05 (2.5 Z ) and the 15 
ages are: 1.0 Ma, 3.16 Ma, 5.01 Ma, 10.0 Ma, 25.12 Ma, 40.0 
Ma, 101.52 Ma, 286.12 Ma, 640.54 Ma, 904.79 Ma, 1.434 Ga, 
2.5 Ga, 5.0 Ga, 11.0 Ga and 13 Ga. These base spectra are dis- 
tributed together with Starlight and reliably fit the stellar popula- 
tions of typical SDSS galaxies ( Cid Fernandes et al.|20 05). The 
age-Z degeneracy is estimated to slightly bias the estimation of 



2010[|Gonzalez Delgado & Cid Fernandes |2010| i. 

We developed a software to adapt the original SDSS spectra 
to the input of Starlight, which performed the following tasks: 
a) uniform re-sampling of the spectra with a final resolution 
of 0.1 nm and transformation from vacuum to air wavelengths 
based on IAU standards (Morton 1991); b) dereddening of the 
flux and its uncertainty based on the dust emission maps of 
Schlegel et al. (1998]) and the extinction curve by O'Donnell 
(1994 , which is an improvement of the extinction curve of 
Cardelli et al. ( 1989|l in the optical and the near infrared; and c) 



shifting of the spectra to the rest-frame using the redshifts from 
the Princeton SDSS data (specBS). The initial values of the ve- 
locity dispersions were also extracted from the Princeton data 
(see Table [2}. 

In Table [3] we present the following physical parameters es- 
timated by Starlight: the velocity shift of the spectrum with re- 
spect to the rest frame, the velocity dispersion of the stellar popu- 
lations, and the extinction. In a few cases (4.7%) Starlight found 
slightly negative extinctions, which may indicate that the deter- 
mination of the underlying stellar populations is not fully reli- 
able; however, the fitted spectra are sufficient to remove the ab- 
sorption features. All fits are clear-cut except for some of the few 
(n < 7) galaxies that present strong broad lines. But the emission 
spectra of these galaxies are not needed for their classification. 
In Table|4]we present the breakdown of stellar populations found 
by Starlight to best fit each spectrum, indicated by the percent- 
age of each different stellar population. Four examples of fits are 
shown in Fig. [4] 

The synthetic spectra obtained from the combination of stel- 
lar populations (as indicated in Table |4]) are available online at 
the AMIGA web page (http://amiga.iaa.es/). 



2.5. Emission lines 

Using the stellar populations derived from the Starlight program, 
we estimated the separate contribution of the emission from the 
interstellar medium, and the AGN emission, if present. The syn- 
thetic spectrum of the stellar populations was subtracted from 
the initial spectrum for each galaxy not flagged in Table|2](galax- 
ies with narrow lines and good data; n = 344), and the narrow 
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Table 2. Catalogue of SDSS nuclear spectra and some of their properties. 1 





plate 


mjd 


fiberld 


Compl. 




z 


V. disp. 
km/s 


z (specBS) 


V. disp. (specBS) Comm. 3 
km/s 


11 


389 


51795 


208 


i 


0.01316 


± 0.00012 




0.013181 ±0.000013 


49.6 ±6.0 


16 


390 


51900 


522 





0.01824 


±0.00015 


103.0 ±2.0 


0.01821 19 ±0.0000054 


102.2 ± 2.0 


19 


753 


52233 


95 





0.01751 


±0.00018 


169.0 ±5.0 


0.0175050 ±0.0000090 


174.6 ±3.8 


56 


400 


51820 


582 


1 


0.017262 


± 0.000099 




0.0172053 ±0.0000044 


134.6 ±4.5 


60 


401 


51788 


223 


1 


0.01727 


±0.00011 




0.017191 ±0.000014 


31 ± 19 


187 


1864 


53313 


171 


1 


0.02767 


±0.00016 


77.0 ±5.0 


0.027633 ±0.000011 


86.0 ±3.9 


189 


1865 


53312 


417 


1 


0.01046 


± 0.00019 


191.0 ±4.0 


0.0104411 ±0.0000067 


189.9 ±3.0 



1 AMIGA galaxies with SDSS DR6 nuclear spectra (n = 362). Columns: (1) CIG catalogue number; (2) SDSS plate; (3) SDSS modified 
Julian date; (4) SDSS fiber Id; (5) code for membership in the complete sample; 2 (6) redshift and its uncertainty; (7) velocity dispersion and 
uncertainty; (8) redshift from the specBS pipeline and uncertainty; (9) velocity dispersion from the specBS pipeline and uncertainty and (10) 
comments. 

2 Code of completeness: 1 if the galaxy belongs to the SDSS complete sample and otherwise. 

3 Comments: "BL" for galaxies with broad lines (see Sec. [5} and "BZ" for galaxies with bad spectral zones (see text, Sect. 2.3 I. 



Table 3. Starlight output. 1 



Table 4. Starlight stellar populations. 1 



CIG 


v shift 


v disp. 


A v 2 


CIG 


Age 


Metallicity 


Percentage 




km/s 


km/s 






Ma 


Z 


% 


11 


-36.82 


55.14 


0.6932 


11 


1 


0.004 


0.0000 


16 


-37.38 


93.08 


-0.0750 


11 


3.16 


0.004 


0.0000 


19 


-37.82 


165.12 


0.1241 


11 


5.01 


0.004 


2.9137 


56 


3.03 


114.42 


0.6307 


11 


10 


0.004 


6.0817 


60 


-16.99 


52.00 


0.5313 










187 


-46.87 


79.07 


0.3914 


11 


11000 


0.050 


1.9735 


189 


-39.79 


183.45 


0.0660 


11 


13000 


0.050 


0.0000 










16 


1. 


0.004 


0.0000 



1 SDSS sample, narrow lines (n = 344). Columns: (1) CIG cata- 
logue number; (2) shift in velocity computed by Starlight; (3) ve- 
locity dispersion of the stellar component; (4) extinction computed 
by Starlight. 

2 In a few cases the computed extinction has slightly negative values 
(see text, Sect. f2A\. 



1 SDSS sample, narrow lines, 45 populations per galaxy (n = 
344 x 45). Columns: (1) CIG catalogue number; (2) age of the stel- 
lar population in megayears (Ma); (3) metallicity; (4) contribution 
of this stellar population. 



lines found in the resulting spectra were measured using a semi- 
automatic Gaussian fit. 

As will be explained in the next section, there are specific 
emission lines that can be used to provide a diagnostic test for 
the presence of an AGN. We therefore performed the fitting in 
five different spectral regions: a) the H/3 region, b) the [Om] 
region, c) the [Oi] region, d) the Ha and [Nn] region, and e) 
the [S n] region. If nearby lines existed (closer than 0.0030 in 
A[log 10 (/l)] units), we used a single fitting region for all nearby 
lines. Those regions correspond to [N n] and Ha and to [S n]. To 
remove any residual continuum emission from the subtraction 
of the underlying stellar spectrum, we subtracted a linear contin- 
uum baseline for each region. An estimate of the best fit to the 
lines in each region was computed by developing an automatic 
pipeline, then the estimated fit was visually checked and changed 
where needed. Line fluxes were calculated using the area below 
the fitted Gaussian: V2~7r (height x width), where width is the 
standard deviation of the fitted Gaussian. The associated error 
is proportional to the noise in the spectral region and is calcu- 
lated as, V27r (cr x width), where cr is the rms continuum noise. 
When the line peak is below 3cr significance, the upper limit is 
calculated as, V2n(3cr x width). The line parameters are pre- 



sented in Table [5] The fluxes of the detected emission lines were 
subsequently compared with those obtained from the MPA-JHU 
SDSS catalogues^] The line flux ratios are consistent within the 
error for all measurements. 



3. Classification and properties of the AGN 

Some galaxies in the SDSS sample show broad emission lines 
and were directly classified as Seyfert 1. Since the Seyfert 
1/Seyfert 2 separation was not the main aim of the work, we 
classified the Seyfert type of the galaxies by visual inspection. 
We found nine Seyfert 1 galaxies (CIG 204, 214, 336, 349, 719, 
747, 749, 893, lOO8)0and among them seven clearly exhibited 
a broad Ha line. CIG 336 and 893 showed a faint broad Ha 
component but, because of the broad width of the [N n] lines we 
classified them as Seyfert 1 . In a future work we will use a para- 
metric classification scheme like the one in |Hao et"aIT1 ( |2005| > and 
|Martfnez etaLlpOOg] ). 



2 http://www.mpa-garching.mpg.de/SDSS/DR7/ 

3 AGN classification obtained from the literature in ISabater et al.l 



( 2008): CIG 214 - Sy 1.0; CIG 349 - Sy 1.5;CIG719-Sy 1.0; CIG 1008 
Sy 1.2 
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Table 5. Measurement of the lines. 1,2 



CIG 


Line 
10- 10 m 


Plux 

l(T 20 Wm- 2 


Pc£ik intensity 
10- 10 Wm- 2 m" 1 


Width 
l(r l0 m 


Pos. 
10- 10 m 


o~ Comment 

10- 10 Wnr 2 nr 1 


11 


Ha (6563) 


566.4 ± 3.7 


136.1487 


1.6597 


6563.50 


0.8954 


11 


[Nil] 6548 


56.5 ± 3.8 


13.2799 


1.6964 


6548.78 


0.8954 


11 


[Nil] 6583 


186.9 ±3.8 


43.9571 


1.6964 


6584.18 


0.8954 


11 


[S ii] 6730 


61.5 ±3.6 


14.7222 


1.6671 


6731.64 


0.8616 


16 


Ha (6563) 


37.3 ±7.6 


5.9584 


2.5002 


6560.72 


1.2170 


16 


[Nil] 6548 


< 22.9 ± 7.6 


1.9993 


2.5002 


6550.32 


1.2170 



1 SDSS sample, narrow lines, eight lines per galaxy (n = 344 x 8). Columns: (1) CIG catalogue number; (2) name of the line. Numerical codes 
in the electronic version: - Ha (6563 x lO" 10 m), 1 - [Nil] 6548 x lO" 10 m, 2 - [Nil] 6583 x lO" 10 m, 3 - H/3 (4861 x lO" 10 m), 4 - [Om] 
5007 x 10~ 10 m, 5 - [Ol] 6364 x lO" 10 m, 6- [Sn] 6716 x lO" 10 m, 7 - [Sn] 6730 x 10" 10 m; (3) flux of the line, the symbol '<' indicates that 
the flux is an upper limit; (4) peak intensity of the line; (5) width of the line; (6) central wavelength of the line; (7) estimated noise for the line 
region; (8) comments about the measurement of the line. 

2 Equivalences of units: 10~ 20 WirT 2 = 10~ 17 erg s~'cnT 2 and 10~ 10 WrrT 2 irT 1 = 10~ 17 erg s~'cirr 2 A -1 . 







Fig. 4. Example of the fit of the underlying stellar populations 
performed by Starlight for four AMIGA galaxies with different 
colours. The original spectrum is shown as a black line, the fitted 
spectrum as a grey line and the line at the bottom indicates the 
unmasked region at which the fit was performed. The average 
ages (log) and metallicities of the estimated stellar populations 
for the depicted galaxies are CIG 657 - 8.8 [a] - 0.005 Z; CIG 
344 - 9.4 [a] - 0.015 Z; CIG 1029 - 9.9 [a] - 0.023 Z; CIG 735 - 
10.1 [a]- 0.020 Z. 

Galaxies presenting only narrow lines were classified using 
the standard emission line diagnostic diagrams (Baldwin et al. 
|1981| |Veilleux & Osterbrock|1987) . Line ratios adopted as suit- 
able diagnostics were a) log([Oni]/HjS), primarily an indica- 
tor of mean ionization level and temperature; b) log([Nn]/Ha), 
which is less immediately obvious but provides a good sep- 



aration between star-forming nuclei and AGN (Osterbrock & 



|Ferland|2006) ; c) log([S n]/Ha) and log([0 i]/Ha), as indicators 
of the relative importance of an extended partially ionised zone 
produced by high-energy photoionization. These ratios were 
grouped to form the three following diagnostic diagrams: 

- log([Om]/H/3) vs. log([Nn]/Ho') (from here-on [Nn] dia- 
gram). 

- log([0 m]/H/3) vs. log([S ii] /Ha) ([S n] diagram). 

- log([Om]/H/?) vs. log([Oi]/Ha) ([Oi] diagram). 



The line ratio measures used in the diagnostic diagrams for the 
AMIGA SDSS sample are presented in Table[6] A ratio is given 
as an upper or a lower limit if only an upper limit was measured 
for one of the lines. When both lines in a ratio involved upper 
limits, the corresponding field is left blank. 



Veilleux & Osterbrock ( 1987 1 proposed several semi 



empirical boundaries within the diagnostic diagrams including 
a theoretical starburst region and objects showing other types 
of excitation such as AGN. LINERs were defined as a sepa- 
rate class of nuclei by Heckman ( 1980 (, and Ho et al. 



1997) 



defined similar classification criteria using the [O i] diagram to 
distinguish between Seyfert 2 and LINER nuclei (see also the 
review of Ho 2008). We used the theoretical curves proposed 
by |Kewley et al.| ( |200"I) to separate narrow line AGN (NLAGN; 
above the curves) from SF nuclei (SFN; below the curves; in the 



[S n] and [O i] diagrams) and the empirical curve of |Kauffmann 
et al. (2003a) to select SF nuclei (SFN; below the curve in the 
[N n] diagram). Galaxies between both curves in the [N n] dia- 
gram were considered to harbour a transition object (TO; com- 
posite of SFN and NLAGN). To separate Seyfert 2 and LINERs, 
we used the empirical classification scheme provided by |Kewley 
et al. ( 2006) for the [Sn] and the [Oi] diagrams. If only one 



of the two line ratios was available for a given diagnostic dia- 
gram, a definitive classification can be obtained from that line 
ratio in some cases ( |Miller et al.|2003[ [Martinez et al.|2008| ), e. 
g. log([Om]/H/?) > 0.8 or log([Nil]/Ha) > 0.0 for AGN. 

Censored data (upper and lower limits) were taken into ac- 
count for the classification. Some galaxies can be unambigu- 
ously classified despite the presence of an upper limit. For exam- 
ple, a galaxy located below the Kauffmann et al. ( 2003a i curve 
(SFN region) is classified as SFN even if there is an upper limit 
in the log([N nJ/Ha) ratio. We can classify 11% more galaxies 
(increasing the overall classification fraction from 84% to 95%) 
if upper limits are considered (see Sect. |3.1| for a more detailed 
comparison). 

The left panel of Fig. [5] shows the [Nn] diagnostic diagram 
for the SDSS complete sample (although classification criteria 
were applied to all galaxies in the SDSS sample). Galaxies clas- 
sified as NLAGN typically have higher errors in their emission 
line ratios because the flux of their emission lines is usually 
lower than for SFN types. In the central and right panels of Fig. [5] 
the [S n] and [O i] diagnostic diagrams are shown. The transition 
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Table 6. Line ratios used for the diagnostic. 1 



7 



CIG 


log([Nn]/Hff) 


log([Om]/H8) 


log([Sli]/Ha) 


log([Ol]/H«) 


11 


-0.4815 ±0.0028 


-0.614 ±0.015 


-0.5814 ± 0.0028 


-1.6132 ±0.0028 


16 


< -0.212 ± 0.089 




< 0.035 ± 0.089 


< -0.053 ± 0.089 


19 


0.025 ± 0.067 


0.07 ±0.13 


0.011 ±0.067 


< -0.217 ±0.067 


56 


-0.31933 ±0.00090 


-0.6795 ± 0.0054 


-0.63006 ± 0.00090 


-1.64935 ±0.00090 


60 


-0.4241 ± 0.0027 


-0.372 ±0.012 


-0.3988 ± 0.0027 


-1.5158 ±0.0027 


187 


0.128 ±0.027 


0.666 ± 0.093 


-0.278 ± 0.027 


< -0.767 ± 0.027 


189 











1 SDSS sample, narrow lines (n = 344). The first column is the CIG catalogue number. The other columns are the different line ratios used in 
the diagnostic diagrams. The presence of an upper or a lower limit is indicated with a '<' or a '>' symbol, respectively. The precision of the 
numbers is adjusted for the error to have two significant digits. 




Fig. 5. Emission line diagnostic diagrams. Error bars have been plotted, and an arrow added for upper or lower limits. Left panel: 
[Nn] diagnostic diagram. The shape of the marker indicates the final classification of nuclear activity: circle - SFN; pentagon - TO; 
triangle - Sy2; square - LINER, diamond - NLAGN. Filled markers correspond to galaxies classified using this diagnostic diagram, 
and the colour indicates the derived classification: red - AGN; green - TO; blue - SFN; cyan - TO or SFN; light green - TO or 
AGN. See Section [3lfor a complet e description of the classification scheme. Middle panel: [S n] diagnostic diagram. The AGN-SFN 
separation curve of |Kewley et al.| ( |2001| ) is shown as a black solid line and the Sy2-LINER separation line of |Kewley et al.| {2006) 
as a green dashed line. Symbol shapes are the same as in the [Nn] diagram and the colours indicate the classification obtained with 
this diagram: red - Sy2; yellow - NLAGN; orange - LINER; and blue - SFN. Right panel: [Oi] diagnostic diagram. Symbol shapes 
are the same as in the [N n] diagram. Colours and separation curves are the same as in the [S n] diagram. Empty symbols correspond 
to galaxies that cannot be unambiguously classified using this diagram but that were eventually classified using the other diagrams 
(see text). 



objects are widely spread in both the NLAGN and star-forming 
regions, as is usually found in the literature. 

To adopt a final classification for each galaxy, we proceeded 
as follows. The classification obtained from the [Nn] diagram 
was given precedence over classifications inferred from other di- 
agrams. In most cases there was agreement (although there are 
nine cases where a galaxy is classified as SFN using [N n] and 
as AGN based on the other diagrams or vice versa). If an ob- 
ject could not be unambiguously classified using this diagram, 
the classification obtained with the other two diagrams was used 
if the two classifications agreed. We could not distinguish be- 
tween TO or NLAGN for three galaxies, and between TO or 
SFN in three other cases. In these cases, we added both classifi- 
cations to the database. Notice that some of the NLAGN cannot 
be sub-classified as LINER or Seyfert 2 using the [S n] and [O i] 
diagrams and were entered as NLAGN in the database. If a tran- 
sition object was classified as an AGN or SFN in the [S n] and 
[O i] diagrams, we indicate this classification in the database in 



an additional "transition object type" column. The final classifi- 
cation is presented in Table [7] in which galaxies that presented 
broad lines have been added as Seyfert 1 . 

The definition of the lines of separation between different 
types of galactic nuclei in the diagnostic diagrams can affect 
the final classification of a galaxy (e.g. see |Kewley et al.||2006| 
and |Constantin & ^V ogeley 2006). Differences in the classifica- 
tion also arise from the method chosen for the final classification 
from the different diagnostics diagrams where these disagree, or 
the different handling of non-detected or upper limits in the fine 
flux. New alternative classification methods (Buttiglione et al. 
|2010[ |Cid Fernandes et al.|2010) > use different combinations of 
lines or derived parameters and classification schemes. Stasinska 
|et al. ( 2008| l studied the possibility that the emission of some 
of the galaxies in the region of AGN could actually be pow- 
ered by stellar processes (hot post-asymptotic giant branch stars 
and white dwarfs; Binette et al.||T994). The classification pro- 



posed by |Cid Fernandes et al.| ( |2010| l allows one to separate be- 
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Table 7. Classification. 1 - 2 



Table 8. Statistics of nuclear activity in AMIGA galaxies. 



CIG 


[Nn] 


[S ii] 


[Oi] 


Final 


TO 




class. 


class. 


class. 


class. 


type 3 


11 


SFN 


SFN 


SFN 


SFN 




16 












19 


NLAGN 


LINER 




LINER 




56 


SFN 


SFN 


SFN 


SFN 




60 


SFN 


SFN 


SFN 


SFN 




187 


NLAGN 


Sy2 




Sy2 




228 


TO 


SFN 


SFN 


TO 


SFN 



SDSS sample (n = 353). Columns: (1) CIG number; (2) classifica- 
tion with the [Nn] diagnostic diagram; (3) classification with the 
[S n] diagnostic diagram; (4) classification with the [O i] diagnostic 
diagram; (5) final classification adopted for our study; (6) transition 
object subtype. 

Acronyms: SFN - star-forming nucleus; TO - transition object; 
NLAGN - narrow line AGN (Seyfert 2 or LINER); LINER - 
LINER; Syl - Seyfert 1; Sy2 - Seyfert 2; UNK - unknown (only 
for the TO subtype). Three galaxies have no unambiguous classifi- 
cation between TO and NLAGN ("TO or NLAGN") and three more 
between TO and SFN ("TO or SFN"), as explained in the text. 
Subtype of transition object if classified as a TO (see text). 



tween "retired galaxies" (powered by these processes) and weak 
accretion-powered AGN ( |Cid Fernandes et al.||2011| l. All these 
differences in the definitions must therefore be taken into ac- 
count to compare between different samples. For this reason, to 
allow a fair unbiased comparison with other samples, we kept 
and release all information (emission line fluxes, detection in- 
formation, emission line ratios, classification for each diagnos- 
tic diagram and final classification) needed to allow a different 
classification scheme to be applied. 

Using our classification criteria, we were able to perform a 
statistical study of the different types of galaxies found in the 
sample. The values obtained are listed in Table [8] both for the 
SDSS sample and the SDSS complete sample, although only the 
latter ones were considered for statistical purposes. The preva- 
lence of optically selected active galactic nuclei for the AMIGA 
sample of isolated galaxies is 20.4%, or 36.7% including TOs. 

3.1. Robustness of the classification method 

The fraction of galaxies that can be classified and their final clas- 
sification may depend on factors such as the signal-to-noise level 
of the spectrum, the classification method used, or the criteria 
used to assign a given type of nuclear activity depending on the 
classification obtained from the diagnostic diagrams. 

We tested the robustness of our classification method with re- 
spect to the signal-to-noise level of the spectra. The required de- 
tection level of the emission lines was raised from 3cr to 5cr, 7<x 
and lOcr (higher <x cut levels mimic a 3cr selection in observa- 
tions taken with progresively lower signal-to-noise ratios). Then, 
we re-calculated the fraction of galaxies of each type using a) 
our method, which takes into account the information carried by 
upper limits and b) a classical classification method, which re- 
quires the detection of the four main emission lines (e.g. Kewley 



et al. 2006). The numbers are presented in Figure|6] The fraction 
of classified galaxies drops with the noise for both methods, but 
the drop is steeper for the classical method and, while with a lOcr 
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19 


7 
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5 


2 
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TO or NLAGN 


1 


1 


0.4% 
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4 
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321 


213 


94.2% 
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191 


130 


57.5% 
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in 

5 1 
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AGN: 


77 


46 


20.4% 


NLAGN total: 


68 


41 


18.1% 


LINER 


27 


18 


8.0% 


Sy2 


10 


9 


4.0% 


NLAGN 


31 


14 


6.2% 


Syl 


9 


5 


2.2% 
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130 


83 


36.7% 
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353 


226 
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Fig. 6. Fraction of galaxies that are a) classified; b) classified as 
SFN; c) classified as TO; d) classified as AGN, with respect to 
the noise level using two different classification methods. The 
first method is the one presented in this paper that takes into 
account the information provided by upper limits (solid line, 
square markers) and the second method is the classical method 
that takes into account only galaxies with detections of all lines 
involved in the [Nn] diagnostic diagram (dashed line, round 
markers). 



cut level our method classifies ~73% of the galaxies, the classi- 
cal method classifies only ~32% of the galaxies. The fraction 
of different types of galaxies remains almost constant with re- 
spect to the noise using our method, with a slight increase of the 
fraction of SFN with high cut levels (lOcr) probably originating 
from the slight decrease of TOs. The classical method tends to 
overestimate the fraction of SFN and underestimate the fraction 
of AGN as the noise increases because only SFN galaxies have 
sufficiently strong emission lines to be detected over the noise. 
These figures suggest that our method not only obtains a higher 
rate of classifications, but also tends to be more robust in the 
classification of different types of galactic nuclei with respect to 
the signal-to-noise level. 
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We also checked the effect of lowering the detection level to 
2cr to include weak emission lines ( |Cid Fernandes et al.|2010| l. 
We risked considering as detections some noise artefacts, but 
we can check if we are excluding many weak active nuclei. The 
overall effect for our sample is a very slight increase in the frac- 
tion of AGN galaxies. The relative fraction of different types of 
AGN over the total number of AGN is similar at 2, 3 and 5cr. 
From 7cr on, there is a drop in the fraction of LINERs (« 70% 
of the classified NLAGN at 2, 3 and 5<x; * 38% at la and 20% 
at lOcr) probably caused by missing retired galaxies with weak 
emission lines (Cid Fernande s et al.|2011| l. 



11.0 1 1 1 1 1 1 1.10 




Morphology 



3.2. Properties of the host galaxies 

The relation between the chance for a galaxy to harbour an AGN 
and its stellar mass has been well-established for different man- 
ifestations of nuclear activity: optical (e .g., |Kauffmann et al. 

2053] ), or X-ray ( |Tasse et al. | 2011) . The 



2003a) , radio ( |Best et al. 



larger the stellar mass, the higher the probability of harbouring 
an AGN. It is also well known that there is a relation between nu- 
clear activity and morphology ( |Moles et al.|1995| ). Accordingly, 
a higher fraction of AGN would be expected in more luminous 
galaxies and in earlier morphological types. 

We studied some properties of the host galaxies with respect 
to the different types of nuclear activity found in the AMIGA 
sample. We considered only galaxies belonging to the SDSS 
complete sample to ensure statistical robustness. The relation of 
L B with L K is almost linear for AMIGA galaxies (Fig. [7J. The 
fraction Lb/Lk has a mean value of 0.946 ± 0.030 with a maxi- 
mum value of 1 .082 and a minimum of 0.884 and presents a clear 
correlation with the morphology, with higher values for late- 
types (Pig. [7j. In Fig.[8]the relation of Lk and the morphological 
type for galaxies in the SDSS complete sample that have NIR 
fluxes available (n — 21 1; ~93.4%) is shown. We show the statis- 
tics for AGN, TOs and SFN as well. Active galactic nuclei are 
hosted by earlier types and higher luminosity galaxies than TOs, 
while just the opposite happens for SFN. Because early types 
present a bigger fraction of the galaxy mass in the spheroidal 
component and AMIGA galaxies harbour pseudo-bulges instead 
of classic bulges (Durb alaet al.|200 9), there might be a relation 
between the prevalence of AGN and the mass of the pseudo- 
bulge in AMIGA isolated galaxies. 

In Fig. [9] we show the fraction of AGN and AGN + TOs 
with respect to optical luminosity (L B ; panel a), L K (proxy for 
stellar mass; panel b) and morphology (panel c). We also con- 
sidered AGN plus TOs plus non-classified galaxies in the case 
of the morphology. There is a monotonic increase of the frac- 
tion of AGN galaxies with respect to Lb and Lk, similar to the 
one found in the literature. The trend is also similar consider- 
ing the transition objects. In terms of morphology, there are no 
AGN or TOs found in very late-type (Scd to Irregular) galax- 
ies. There is a monotonic increase in the fraction of AGN (from 
0% to 100%) from type Sc to Sa, while earlier-types yield nois- 
ier results because there are few galaxies in these bins. There 
are four early-type galaxies presenting an SFN. Three of them 
(CIG 364, 498 and 620) are edge-on galaxies classified as SO or 
SO/Sa, and CIG 503 is an SO galaxy presenting Wolf-Rayet sig- 
natures ( |Brinchmann et al.|2008) a nd identified as an early-type 
star-forming galaxy by Temi et al. ( 2009) and Wei et al. ( 2010) . 
The increase of the fraction of AGN is steeper with respect to 
Lk than with respect to Lb due to the correlation between Lb/Lk 
and the morphological type: galaxies with higher L K and later 
type tend to have lower Lb and vice versa, thus, reducing the 
steepness of the slope (see Fig. [7). 



Fig. 7. Luminosities of the AMIGA isolated galaxies that are 
in the SDSS complete sample and have available NIR lumi- 
nosities (n = 211). The left panel shows the relation between 
Lb and Lk- The ordinary least-squares (OLS) regression line is 
shown as a solid black line. The right panel shows the relation of 
log(LB)/ log(Lic) with the morphology. The OLS regression for 
the whole sample is plotted as a solid line. Two additional OLS 
regression lines are shown as dashed lines, the first for galaxies 
with morphological types between E and Sbc and the second for 
the rest. A small random scatter (from -0.25 to 0.25 of a class) 
was added to the morphology values to allow a proper visualisa- 
tion of the symbols. 



Seyfert 1 galaxies tend to present higher luminosities than 
narrow-line AGN galaxies, as shown in Fig. [8] This difference 
may be produced by weak broad lines that are missed in galax- 
ies with lower luminosities because of the lower signal-to-noise 
level or the direct detection of the AGN light in the luminos- 
ity measurements. There is no clear dependence on luminosity 
for the other types of AGN, but the fraction of LINERs relative 
to the number of classified NLAGN (LINER + Sy2) depends 
strongly on the morphological type. It ranges from 25% for Sc 
types to 83% for Sab and earlier types (Sc - 25%; Sbc - 50%; 
Sb - 82%; Sab and earlier - 83%). This trend is consistent with 
some of the LINERs being actually retired galaxies ionized by 
low-mass evolved stars (Cid Fernandes et al. 2011). 



4. Comparison with galaxies in isolated denser 
environments 

One of the main aims of this study is to determine the effect of 
the environment and interactions on the prevalence of nuclear 
activity. After estimating the normal level of nuclear activity for 
a sample of galaxies not affected by major interactions, a care- 
ful comparison with samples of galaxies in denser environments 
is required to quantify the effect of the environment/interaction. 
An unbiased comparison will require to take into account the 
effect of the density-morphology and the density-luminosity re- 
lations (|Dressler|[T980l [Kauffma nn et al.||2004[ |Blanton et al.| 
2005 1 |Deng et al.||2009| |Park & Choi||2009| |Deng et al.||2011) . 
The strong correlation of the fraction of AGN with morphology 
and luminosity (both optical and NIR) could otherwise produce 
a strong bias when comparing with samples that could be com- 
posed of galaxies with different morphologies and luminosities. 

Some of the densest environments in the universe are in com- 
pact groups of galaxies, with projected galaxy densities simi- 
lar to those found in the cores of dense clusters. Among them, 
Hickson Compact Groups (HCG; characterized by 4 to ~ 1 



members with a low velocity dispersion ~ 200 km s 
1982) have been selected using an isolation criterion 



Hickson 



Hickson 
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Fig. 8. Properties of the AMIGA isolated galaxies that are in the 
SDSS complete sample and have available NIR luminosities (n 
= 211). The main panel shows L K versus the morphology. The 
SFN are marked as white round symbols, unclassified galaxies 
as black + symbols, TOs as cyan squares, and AGN as trian- 
gles (red triangles in general and yellow inverted triangles for 
the Seyfert 1). A small random scatter, like in Fig. [7] was added 
to the morphology values. The additional panels show the statis- 
tics for different types of nuclear activity with respect to the mor- 
phology (upper panel) and Lk (right panel). The box marks the 
lower and upper quartiles, the median is indicated as a red line 
inside the box, the whiskers extend to the most extreme datapoint 
within 1 .5 x (xj5% - ^25%) from the median, and the outliers are 
marked as plus symbols. The mean values are marked using a 
star symbol. The mean of L K for Seyfert 1 galaxies is indicated 
in the right panel with a cross at the left of the AGN box. 



[19821 |Sulenfic|1987||Hickson et al.|1992"l >. This combination of 
high galaxy density in a low-density environment makes them 
a unique sample for comparison with isolated galaxies by min- 
imising the possible effects of the large-scale environment. It is 
important to remark that AMIGA and HCG samples were de- 
fined using only isolation criteria. 

Martinez et al. ( |2010 1 studied the optical nuclear activity of 



a complete well-defined sample of HCG galaxies (n = 270) ob- 
taining new spectroscopic measurements (n = 200), as well as 
archival spectra (n = 11) and emission line measurements from 
the literature (n = 59). The spectra were corrected for absorp- 
tion features from the underlying stellar populations, further- 
more, the authors followed a classification scheme very simi- 
lar to the one that we present in this paper: the diagnostic di- 
agram used is the [Nn] one, separation lines are identical to 
ours and the log([N 11] /Ha) line ratio is used when it is the only 
one with confident values. The mean signal-to-noise ratios of 
the log([N ii]/Hq-) coefficient for both samples are very similar 
(HCG - ~ 42; AMIGA - ~ 56). 

The origin of the luminosity and morphology data used by 
Martinez etaLl ( |2008[ ) is HyperLeda ( |Patureletal.|2003| ). We also 



Fig. 9. Relation between the fraction of AGN (red solid line, 
square symbols) or AGN plus TOs (green dashed line, round 
symbols) with the optical luminosity (a), NIR luminosity (b), 
and morphology (c) for the SDSS complete sample. In the mor- 
phology panel we also show the fraction (black dotted line, di- 
amond symbols) of AGN plus TOs plus unclassified galaxies 
(galaxies without detected emission lines that are often referred 
to in the literature as passive galaxies). 



catalogue for AMIGA galaxies. AMIGA curated data (presented 
in Sect. |2.2| > present some fundamental differences compared to 
Hyperleda data, which render a direct comparison between them 
unsuitable. In particular, the optical luminosity is corrected for 
different effects (Fernandez Lorenzo et al. 2012 1 and presents 
a systematic shift with respect to Hyperleda data. The differ- 
ence in morphological classification, although very low on av- 
erage (-0.6), presents a high dispersion (2.5). Hence, we used 
HyperLeda data for the comparison between AMIGA and HCG 
to allow an homogeneous comparison. The morphological clas- 
sification is coded into numbers using the RC3 morphologi- 
cal types system and the absolute magnitude was transformed 
into optical luminosity (B-band) using the relation log(LB) = 
(5.51 - M B )/2.5[L©]. Finally, there are data for all but two 
AMIGA galaxies of the SDSS complete sample (n = 224) and 
for all but one galaxy of the HCG sample (n = 269). The ranges 
and means of the optical luminosity are very similar for both 
samples (AMIGA: 10.31 ± 0.49; HCG: 10.4 ± 2.0) but the dis- 
tribution of morphologies is quite different (AMIGA: 3.6 + 2.5; 
HCG: 0.4 + 3.6) with earlier types in compact groups (see upper 
panels of Fig. 10 1. We will present a method to take this different 
distribution into account. 



obtained luminosity and morphology data from the HyperLeda 



4. 1 . Comparison method 

An empirical probability density function (pdf) is built that con- 
siders the optical luminosity, morphology, and the classification 
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as an AGN for a sample of galaxies. The definition is as follows: 



f(L,t,a) = - V 



1 



2no- L icr t i 



exp 



2 <i 



2\ 



exp 



(f ~ tj) 

2 <- 



2 A 



where L is the luminosity, f the morphology in the RC3 system, 
and a is a value that indicates whether the galaxy harbours an 
AGN (a = 1) or not (a =0). Note that in the summation term, a 
galaxy i is only added if a and a, are equal. 

The partial probability density function with respect to the 
AGN activity represents the distribution in luminosity and mor- 
phology of the sample and is defined as f(L, f) = £ a f(L, t, a) = 
f(L,t, 1) + f(L,t,Q). We can obtain for each L and t pair the 
mathematical expectation with respect to the values of a: 

E[a](L, t) = f(L, t, l)/f(L, t) = p a (L, ?), 

Hence, the probability for a galaxy to harbour an AGN for a 
given luminosity Lq and morphology to, derived from our em- 
pirical pdf, is: p a= i(Lo,to). To compare the prevalence of AGN 
in two different samples, we can compare the actual number of 
AGN galaxies in the reference sample with the estimated num- 
ber of AGN galaxies for a second sample calculated using the 
same p a= \(L, t) but weighted by the actual distribution of L and 
t in the second sample. In other words, the estimated number of 
AGN galaxies in a region (R) of the L, t plane for the second 
sample (s2), using the distribution of the first sample (si), is 



n s 2,AGN 



ns2 K f > 



s2 (L,t)p sl . a=l (L,t)dLdt. 



In that way, we can check whether estimated and actual values 
are compatible within the error (or not), in a given region, taking 
into account the effect of the luminosity and morphology distri- 
bution. 



4.2. Comparison with galaxies in compact groups 

There are 46 (20.5%) galaxies classified as AGN in the AMIGA 
sample and 66 (24.5%) galaxies classified as AGN in the HCG 
sample. The actual distribution of galaxies and AGN are shown 
in the upper panels of Fig. 10 Clearly, the distribution of galax- 



ies is different in the two samples with later types in the AMIGA 
sample and earlier types in the HCG one. The distribution of 
AGN also looks different, but this may be caused just by the 
general distribution of galaxies in the sample. Hence, we ap- 
plied the method described above to discard this possible bias. 
The number of AGN was compared in the region in common 
between the two samples defined by the intersection of regions 
in which 95.45% of the galaxies are within the area for each 
sample (equivalent to 2 sigma in a normal distribution). This 
common region contains « 206.8 AMIGA galaxies and m 221.2 
HCG galaxies and is shown as a contour on the panels of Fig.fTO] 
Values will be measured within this common region from now 
on. The fraction of Seyfert galaxies with respect to the total num- 
ber of AGN in HCG is 37% ( [Martinez et al.|2010| >, while we find 
a fraction of ~ 30% in AMIGA galaxies. Therefore, we do not 
expect to be comparing very different types of AGN (or retired 
galaxies) between the two samples. 

Table [9] shows the number of AGN in each sample. It also 
shows the number of AGN that would be obtained by combining 
the galaxy distribution of each sample with the AGN probabil- 
ity distribution of the other sample (as defined in Sect. 4.1 1. A 
binomial test comparing the observed with the cross-estimated 
number of AGN galaxies yields p-values well above the limit of 




o 


9 




4 




12 ■ 




11 


hi' 




13 


10- 


d. 


to 
o 


9 





















c) AMIGA AGN 



— 11 




5m E 



SO/a Sc 
Morphology 



Sm E 




SO/a Sc 
Morphology 



Sm 



Fig. 10. Distribution of galaxies (morphology and optical lumi- 
nosity) and comparison of the actual and estimated distribution 
of AGN in the AMIGA and HCG samples. In the top panels the 
total distribution of galaxies for the AMIGA (a) and HCG (b) 
samples are shown. The colour ranges from red to yellow, with 
red corresponding to a higher density of galaxies. In the middle 
panels the distributions of AGN in the AMIGA (c) and HCG (d) 
are shown. In the lower panels, the cross-estimated distributions 
of AGN in AMIGA (e) and HCG (f) are shown, as produced 
by combining the actual distribution of galaxies with the proba- 
bility of harbouring an AGN determined from the other sample. 
The colour scale of the middle and lower panels is one third of 
the scale of the corresponding upper panel. The black contour 
marks the region in common between both samples (95.45% of 
the galaxies inside the contour for each sample) and the rectan- 
gles shown in the panel "c" denote the bins chosen for the x 1 
test. 



a — 0.1 for both samples. Hence, we cannot reject the hypoth- 
esis that the AGN fractions are the same, after correcting for 
the luminosity and morphology distribution. The lower panels 
of Fig. 10 show the cross-estimated distribution of AGN apply- 
ing our method. The cross-estimated distributions (e and f) are 
not very different from the observed ones (c and d respectively), 
as is confirmed by the statistical tests. There is a slight discrep- 
ancy in the distribution of AGN for early spirals between the 
observed and estimated distribution in HCG (panels d and f) but 
the low number of galaxies in this region makes this difference 
less significant. 

We also compared the prevalence of AGN in both samples 



using a x test- We were limited by the relatively low number 
of AGN in both samples, therefore, the bins were selected to 
allow a significant number (n > 6) of observed and expected 
AGN in each one: a) Morphology from E to Sa and log(L B ) 
from 9.5 to 11.0, b) morphology from Sab to Sd and log (Lb) 
from 9.5 to 11.0, and c) the region out of the other two bins. 
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Table 9. Comparison between AMIGA and HCG samples. 



Number of objects 
p AMIGA p HCG 



p val." 



AGN 



AMIGA 44.2 
(Fig. 

HCG 56.5 
(Fig. 



10;) 
76.2 



AGN + TO 



SFN 



AMIGA 
HCG 
AMIGA 
HCG 



10) 



79.3 
95.5 - 

71.6 : 



T9~ 

7.2 
VTA 

5.7 



52.2 
(Fig. 
56.0 
(Fig. 



: 6.0 



10;) 
75.4 



52.2% 
46.0% 



101) 



108.2 : 

100.9: 

98.1 ± 
43.6 h 



TT 
7A 

5.2 



5.9% 
16.7% 
0.05% 
6.9 x 10~ 7 % 



1 The rows correspond to the distribution of galaxies (sample) and the 
columns to the probability of harbouring an AGN, an AGN + TO 
or a SFN. The numbers in italics are the cross-estimated, calculated 
using the presented method. 

2 The p-value is computed using a binomial test for a given probabil- 
ity and two different distributions. 



The bins are shown in the left-middle panel of Fig. [TU] The p- 
values for the comparisons are 0.181 and 0.165, both above the 
0.1 limit, hence, AGN distributions in AMIGA and HCG are 
indistinguishable. 

If we take into account both AGN and TOs, the p-value de- 
creases (see Table |9]l and there is increasing evidence of a dif- 
ference between the distributions. This difference could be ex- 
plained by the higher weight of the star formation component in 
TOs. 

Finally, the fraction of SFN nuclei in both samples was 
compared. There is evidence of a low star formation rate 



( Verde s-Montenegro et al. 1998, Iglesias-Paramo & Vflchez 



999 ,|Coziol et al.|2000)> or even suppression of star formation 
( Sulentic et al.||2001 Durbala et al. 2008a) in compact groups. 
We checked if the rate of star formation was different after cor- 
recting for the morphology luminosity distribution. The result 
is shown in Table [9] In this case, we can reject the hypothesis 
that the prevalence of SFN in both samples is similar. We obtain 
always lower values in HCG even after the correction. 



5. Discussion and conclusions 

We presented a careful estimation of the optical nuclear activ- 
ity in a well-defined, statistically significant sample of isolated 
galaxies using SDSS data. We aimed to determine the prevalence 
and properties of nuclear activity in galaxies not affected by ma- 
jor tidal interactions during the last part of their lifetime. To shed 
light on the effect of interaction on nuclear activity, we compared 
our sample with a sample of galaxies that are strongly affected 
by interaction (HCG galaxies). 

From the 6th Data Release of the Sloan Digital Sky Survey, 
we obtained spectral data that were inspected in a semi- 
automatic way to ensure the accuracy of the data. We subtracted 
from the spectra the underlying stellar populations obtained us- 
ing the software Starlight and fitted the emission lines. We ap- 
plied the typical diagnostics diagrams to classify the type of nu- 
clear emission using a classification scheme that takes into ac- 
count censored data. A catalogue of spectroscopic data, stellar 
populations, emission lines and optical nuclear activity classifi- 
cation was provided. The entire intermediate data are released to 
allow the use of different nuclear activity classification criteria in 
the future. Finally, we presented a method to compare the preva- 
lence of AGN with galaxies in denser environments that takes 



into account the effect of the density-morphology/luminosity re- 
lations to avoid a bias in the comparison. The prevalence of 
AGN and SFN was compared with galaxies in Hickson Compact 
Groups. 

From the present study we found that 

1. Our nuclear activity classification method, which takes into 
account the information carried by censored data, seems to 
be robust with respect to an increase of the noise level. 

2. The prevalence of optical nuclear activity in AMIGA iso- 
lated galaxies is 20.4%. This percentage considers Seyfert 
2, Seyfert 1, LINERs, and unclassified NLAGN (LINER or 
Seyfert 2). If we consider these types together with TOs, the 
percentage rises to 36.7%. 

3. The fraction of AGN increases steeply towards higher Lk 
(proxy for stellar mass) or Lb and earlier morphological 
types. 

4. We found no evidence of a difference between the pop- 
ulation of AGN in AMIGA galaxies and in Hickson 
Compact Groups after correcting for the effect of the density- 
morphology/luminosity. On the other hand, we found signif- 
icant evidence of a difference for SFN with a higher preva- 
lence in isolated galaxies with respect to galaxies in HCGs, 
which is consistent with a lowered SF in HCG. 

The isolation criteria for AMIGA galaxies imply nearest- 
neighbour crossing times of at least 3 Ga (|Verdes-Montenegro| 



etaLl[2003trVerley et al.|2007al [Argudo-Fernandez et al.|2012 



in prep.) That fact, and the general properties of the sample, leac 
us to consider that AMIGA galaxies were not affected by ma- 
jor tidal interactions during the last part of their lifetime. The 
time-scale of the duty cycle of an AGN is very short (less than 
100 Ma; |Haehnelt et al.|[T998l |Greene & Hol|2007t [Ho][2008 ; 



Schawinski et al.|2009| > in comparison with the isolation time of 



our galaxies. Hence, AGN activity in AMIGA galaxies should 
be unrelated to a major interaction event. 

Major interactions have been proposed as a mechanism to 
trigger AGN activity in some cases (se e review by |Combes[2 003 ; 
Ellis on et al.||20lT Liu et al. 2012 1, although no clear rela- 



tion was found in some studies using the SDSS ( Li et al.j 2006 



2008 ). Seyferts do not require external interactions to be fuelec 
(e.g. |Comb es 2003), thus, alternative mechanisms such as minor 
mergers, accretion of gas from cosmic filaments and/or internal 
mechanisms (Combes 2006) are needed to explain the presence 
of optical active nuclei in isolated galaxies. On the other hand, 
these alternative mechanisms do not trigger radio nuclear activ- 
ity in isolated galaxies and an additional factor, linked to the en- 
vironment (directly or indirectly), is needed to explain the pres- 
ence of a radio jet ( |Sabater et aLl|2008[ |Sabater||20"09"t |Sabater 
|etal.|2010p . 

If we consider a traditional feeding mechanism of black 
holes driven by major mergers (Korm endy et al.|2011[ ), the rel- 
atively high fraction of AGN in isolated galaxies, dominated by 
pseudo-bulges (Durbala et al. 2009), might favour an alterna- 
tive local feeding mechanism driven by the slow accretion of 
cold gas (Hopkins & Hernquist 2006 ). This mechanism is prob- 
ably closer related to the activity found in galaxies with pseudo- 
bulges piang et al.|201 1[ ). This could mean that the black holes 
found in AMIGA galaxies are close to the black hole seeds pro- 
posed by Kormendy et al. (2011 (, the seeds of the larger black 



holes that are grown by mergers. 

The comparison with galaxies in compact groups yielded 
that the difference in the fraction of AGN between isolated 
galaxies and those in compact groups is not significant (after cor- 
recting for the morphology and luminosity distribution), while 
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for SFN it is clearly significant. However, we recall that we are 
comparing fractions and not intensities or types of SF or AGN 
activity. One possible explanation is that the expected increase of 
AGN activi ty in HCG due to interactions (El lison et al.|201 l||Liu| 



environmental density ( 


Kauffmann et al.|2004| von der Linden 


etal.|2010||Gavazzi eta 


.|201 1). An additional consideration that 



has to be taken into account is the possible contamination of the 
different samples by retired galaxies, whose nuclear emission re- 
sembles the emission of weak AGN (usually LINERs) but is ac- 



tually powered by low-mass evolved stars (Cid Fernandes et al. 
2010 201 1) . To clarify those points, the results of this paper will 



be extended in the future to study the independent effects of en- 
vironment and interaction on the triggering of both optical and 
radio nuclear activity in a large homogeneously selected sample 
of SDSS galaxies. 

Even though isolated galaxies do not present radio nuclear 
activity, the prevalence of optical nuclear activity is not negli- 
gible. Indeed, this leads us to conclude that a major interaction 
is not a necessary condition for the triggering of optical nuclear 
activity. Hence, at least part of the growth of a black hole could 
be produced without the need of a major merging event. 
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Table A.2. Hyperleda data for AMIGA galaxies. 1 



AMIGA 



AMIGA 



1050 galaxies 



Table 1 



5DSS DR6 photometric 
coverage 



AMIGA SDSS phot. 



549 galaxies 



Table A.l 



5DSS DR6 spectroscopic 
coverage 



AMIGA SDSS spec, 
(raw) 



362 galaxies 



Table 2 



corrupted spectra 
removed 



AMIGA SDSS spec. 
Total sample 



353 galaxies 



Table 7 




AMIGA SDSS spec, 
narrow lines 



344 galaxies 



Tables 3, 4, 5 and 6 



complete sample 
selected 

AMIGA SDSS spec. 
Complete sample 



226 galaxies 



Table A. 2 



HCG 



270 galaxies 



Table A. 3 



Fig. A.l. Description of the samples and data presented in the 
different tables within the paper. 



Appendix A: Additional data 

There are some additional data associated with this study: a) 
SDSS photometric data that were processed but not used for this 
study and b) Hyperleda data that were compiled for the com- 
parison presented in Sect. [4] These data are presented here and 
in http://amiga.iaa.es/ for future reference and to allow the data 
replication and reproducibility. 

A diagram of the samples and data presented in the paper is 
shown in Fig. 



A.l 



A.1. SDSS photometric data 



We cross-correlated the positions of our galaxies ( |Leon &] 
|Verdes-Montenegro|2003| l with the SDSS DR6 photometric cat- 
alogue using the interface available on the SDSS web-page 
(search radius of 0.5')- We obtained a list of photometric sources 
near the estimated centre of each galaxy as well as the corre- 
sponding SDSS spectral data when available. The right objects 
were identified using a semi-automatic selection algorithm se- 
lecting the object with the highest z band flux followed by a vi- 
sual inspection as a cross-check. Only ~ 8% of the photometric 
objects were clearly misclassified with the usual source of con- 



CIG 


Morpho. 


M B 


11 


5.1 


-20.132 


56 


3.4 


-20.641 


60 


2.8 


-19.503 


187 


3.9 


-20.805 


189 


-5.0 


-19.511 



1 SDSS complete sample (n = 226). Columns: (1) CIG catalogue 
number; (2) morphological classification (RC3); (3) absolute mag- 
nitude in B-band. 



Table A.3. Hyperleda data for HCG galaxies. 1 



HCG 


Morpho. 


M B 


class. 


HOla 


4.6 


-21.768 


SFN 


HOlb 


-2.9 


-20.564 




HOlc 


-2.3 


-20.214 


AGN 


HOld 


-1.5 


-19.155 


AGN 


H03a 


4.3 


-20.635 


AGN 



1 HCG galaxies (n = 270). Columns: (1) HCG catalogue name; 
(2) morphological classification (RC3); (3) absolute magnitude in 
B-band; (4) nuclear activity classification. Possible values: AGN, 
SFN, TO, AGN/TO or unclassified (-). 



fusion involving a nearby star. A total of 549 AMIGA galaxies 
were included in DR6. 

In six cases the photometry for the galaxy was assigned to 
an object whose estimated position was shifted from the cen- 
tre of the galaxy. In these cases, we manually chose the correct 
photometric object. We inspected in more detail some special 
cases that were finally discarded from our catalogue (CIG 402 is 
strongly affected by the light of a nearby star and CIG 802 and 
388 are resolved Milky Way satellites). As a final check we com- 
pared the apparent magnitudes in g-band from the SDSS data 
and in B-band from the AMIGA database ( |Verdes-Monten egro 
|etal.|2005] l. Although the filter shapes are different and the shape 
of the spectral energy distribution of the galaxy directly affects 
their ratio, a high difference between these two magnitudes may 
indicate an error in the selection of the photometric data. The 
mean of the ratio is 1 .044 and the standard deviation is 0.070. In 
seven cases the difference was found to be above 3cr, suggest- 
ing that the SDSS photometric data might refer to a restricted 
region rather than the entire galaxy. The photometric data for 
these cases were flagged. The values for the SDSS photometry 
are presented in Table |A.l] 

A.2. Hyperleda data 

Hyperleda data were compiled for the AMIGA and HCG sam- 
ples (see Sect. |4j. In Table A. 2 the compiled morphological 
classification and absolute magnitude for AMIGA galaxies are 
shown. In Table |A.3| the compiled morphological classification 
and absolute magnitude plus the nuclear activity classification 
obtained from Mar tinez et al.| ( [20T0) l are shown. 



J. Sabater et al.: The AMIGA sample of isolated galaxies, Online Material p 2 
Table A.l. Catalogue of SDSS photometry. 1 



CIG 


Objld 


u 


g 


r 


i 


z comment 


11 


588015508197212201 


15.76 


14.45 


13.83 


13.52 


13.33 


12 


588290879491408031 


17.04 


15.57 


14.88 


14.55 


14.30 


16 


587731187816857647 


17.04 


15.26 


14.51 


14.12 


13.83 


19 


587730773354479638 


16.83 


14.95 


14.10 


13.68 


13.32 


33 


587731186208669747 


14.90 


13.57 


12.95 


12.65 


12.40 



1 AMIGA galaxies with SDSS DR6 photometry (n = 549). Columns: (1) CIG catalogue number; (2) SDSS Objld of the selected SDSS 
photometric object. In some cases it does not match the Objld associated to the spectroscopic object. (3) to (7) SDSS photometric magnitudes; 
(8) comments for galaxies with poor photometry (see text). 



